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Abstract. In this paper we consider two inverse problems on a closed con- 
nected Riemannian manifold (M, g) . The first one is a direct analog of the 
Gel'fand inverse boundary spectral problem. To formulate it, assume that M 
is divided by a hypersurface £ into two components and we know the eigenval- 
ues Xj of the Laplace operator on (M, g) and also the Cauchy data, on S, of the 
corresponding eigenfunctions cf>j, i.e. (f>j\s, d v (f>j\s, where v is the normal to £. 
We prove that these data determine (M, g) uniquely, i.e. up to an isometry. 
In the second problem we are given much less data, namely, Xj and 4>j\T, on ly- 
However, if £ consists of at least two components, Ex, £2, we are still able to 
determine (M, g) assuming some conditions on M and S. These conditions are 
formulated in terms of the spectra of the manifolds with boundary obtained 
by cutting M along Si, i = 1, 2, and are of a generic nature. We consider also 
some other inverse problems on M related to the above with data which is 
easier to obtain from measurements than the spectral data described. 



1. Introduction and main results 



In this paper we consider some inverse spectral problems on a closed connected 
Riemannian manifold (M,g). The first motivation to consider inverse problems 
on Riemannian manifolds comes from spectral geometry. The famous problem 
here, posed by Bochner and formulated by Kac in the paper "Can one hear 
the shape of a drum?", [14], is the problem of identifiability of the shape of 
a 2— dimensional domain from the eigenvalues of its Dirichlet Laplacian. More 
generally, the question is to find the relations between the spectrum of a Rie- 
mannian manifold (M,g), i.e. the spectrum of the Laplace-Beltrami operator 
— A g on it, and geometry of this manifold. In particular, one can ask, following 
Bochner-Kac, if the spectrum of — A g determines the geometry. However, already 
in 1966, it was known that, in higher dimensions, the answer to this question is 
negative. Indeed, in 1964 Milnor [27] found the first counter-example, a pair of 
Riemannian flat tori of dimension 16 which are isospectral, i.e. the spectra of 
their Laplacians coincide, but not isometric. As for the original Bochner-Kac 
problem in dimension 2, the answer was found only in early 90th. Namely, in 
1985 Sunada [32] introduced a method of producing examples of non-isometric 
isospectral compact connected Riemannian manifolds. Although in this paper 
Sunada did not give the answer to the Bochner-Kac problem, in 1992 Gordon, 
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Webb and Wolpert [10] extended Sunada's method and settled in the negative 
this famous problem by constructing two simply connected non-isometric plane 
domains which are isospectral. Since then there was much activity in this direc- 
tion with many beautiful mathematical results regarding relations between the 
spectrum and geometry, see e.g. [11] and [36] for the current state of art in this 
area. 

It is clear from the above that, in order to determine geometry of a closed Rie- 
mannian manifold, further spectral information is needed. The nature of this 
information can be found if we look at inverse boundary problems. In this case, 
the data given is the trace on dM of the resolvent of the Laplacian with, say, Neu- 
mann boundary condition. Depending on whether this information is given for 
one or many values of the spectral parameter, these inverse boundary problems 
were originally posed by Calderon [7] and Gel'fand [9]. These inverse boundary 
problems were solved, at least on the level of uniqueness and sometimes condi- 
tional stability, for the Laplace-Beltrami operator and also the other types of 
scalar operators, in e.g. [l]-[6], [16, 18, 19], [23]-[25], [28]-[30], [33, 34] and 
monographs [13] or [17] with further references therein. 

As the first inverse problem considered in this paper is an analog of the Gel'fand 
inverse boundary problem, we first reformulate the Gel'fand problem in an equiv- 
alent form which, however, has more "spectral" flavor. Namely, let Aj and <pj be 
the eigenvalues and normalized eigenfunctions of the Laplace operator with Neu- 
mann boundary condition, 

(-A g - = in M, d^jldM = 0; {4>j,4>k)L 2 (M) = S jk , (1.1) 

where d u is the normal derivative to dM. Then the Gel'fand problem [9] is 
the one of the determination of (M, g) from the boundary spectral data, i.e. 
dM, {Xj, <f>j\dM}jLv Note that, due to the Neumann boundary condition in (1.1), 
we do actually know the whole Cauchy data <j)j\dM, d u (f>j\gM of the eigenfunctions 
on dM. To formulate its analog for a closed connected manifold, assume as ear- 
lier that we deal with the Laplace operator — A g on, now closed, Riemannian 
manifold (M,g). As earlier, we assume that our spectral data are given on an 
(m — 1)— dimensional submanifold £ C M, £ = dS, where S C M is an open 
set. We note that this type of data is natural for various physical applications 
when sources and receivers are located over some surface in space rather then are 
scattered over an m— dimensional region or put on, probably remote, boundary 
of M. Such localization is used e.g. in radars, sonars, and in medical ultrasound 
imaging when a single antenna array is used to produce the wave and to measure 
the scattered wave. It is typical also in geosciences/seismology where sources and 
receivers are often located over the surface of the Earth or an ocean. 
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Returning to the possible generalizations of the Gel'fand problem to inverse prob- 
lems on closed manifolds and taking into account that £ splits M into two man- 
ifolds with boundary, S and M \ S, we believe that the most straightforward 
generalization of the Gel'fand problem is the following: 

Assume that we are given the Cauchy spectral data, 

{E, (Aj,0j|e,9i/^-|e)^=i} 

Does these data determine (M,g)1 
Then the first main result of the paper is 

Theorem 1.1. The Cauchy spectral data (1.2) determine the manifold (M,g) up 
to an isometry. 

Probably a more surprising result concerns with another inverse spectral problem 
associated with £. In this case we have less data, namely, the Dirichlet spectral 
data, 

{£, (Xj, 0j|e)£1i}, (1-3) 

and would like to determine (M, g) from these data. It turns out that this is still 
possible when the set S consists of two subsets, 

s = s 1 us 2 , 5!n5 2 = 0, 

if we assume some generic conditions on S\,S 2 . To state this result, we define 
the Dirichlet spectrum of the Laplace operator on a manifold S with boundary 
to be a collection of all eigenvalues Aelof — A g such that 

-A g (f) = X(f), in S, (f)\ d s = 

with some non-zero 0. 

Theorem 1.2. Assume that the Dirichlet spectra of the Laplace operators on S±, 
S 2 , M \ S, M \ Si and M \ S 2 are all disjoint. Then the Dirichlet spectral data 
determine the manifold up to an isometry. 

The paper is organized as follows. In Section 2 we present some auxiliary results 
for the transmission problems on M and the metric reconstruction on E. Section 
3 deals with the reconstruction of the Riemannian manifold (M, g) from the 
Cauchy spectral data and contains the proof of Theorem 1.1. Section 4 is devoted 
to the inverse problem with the Dirichlet spectral data and provides the proof of 
Theorem 1.2. Section 5 contains some further generalizations of the problem and 
its alternative formulations which may be useful in practical applications. 



(1.2) 
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2. Auxiliary results 



2.1. Transmission problem at fixed frequency. Consider a closed connected 
smooth m— dimensional Riemannian manifold (M, g) and an open non-empty 
set S C M with smooth (m — 1)— dimensional boundary Y, := dS ^ 0. Let 
u : M — > K. be a function such that its restrictions w + , w„ onto M + — M\S and 
M_ = 5 are iJ 2 -smooth. We define the traces of u on the different sides of E by 

(w±|s) (x) — lim u(x ± hv(x)), i6E, (2.1) 

(ctw±|s) (x) = lim d v u{x ± hu(x), x e S, 

where z/(x) is the unit normal to E at rr pointing towards M_ and the limits are 
understood in the sense of traces in Sobolev spaces. Denote by [u] = u + \^ — u_\y, 
and [d u u] = d v u + \^ — d v U-\z the jumps of u and its normal derivative across E. 

Consider the transmission problem 

(-A g - \)u := -g-^ditfl^dju) - \u = in M \ E, 
[u] = / on E, [d u u] — h on E, 

where (7 = det((?jj) and [g*- 7 ] is the inverse matrix of [gij]- Although, in general, 
problem (2.2) may be considered with / e # 3 / 2 (E), h € # 1/2 (E), it will be 
sufficient and convenient for us to take f,h<E C°°(E). Problem (2.2) is equivalent 
to the problem 

(-A, - A)« = hfa + /cW s , (2.3) 
where h5s and fd v 5-£ are distributions defined as 



(M S ,V) = / W\*dS g , (fd v 8^) = - / /(cW) |s^, 
dS'p being the volume element on E, for any ip e C°°(M). 

Let A ^ a(—A g ), where a(—A g ) is the spectrum of the Laplace operator — A g . 
Then problem (2.2) has a unique solution. Moreover, its formulation in form 
(2.3) makes it possible to represent this solution, u = u^ ,h (x, A), for x G M \ E, 
as a sum of a single- and double-layer potentials, 

u f ' h (x, A) = u°' h (x, A) + u f '°(x, A), (2.4) 



u 



0,h 



(x,X) = S x h, u f '°(x,X) = -D x f, 



S x h(x)= / G(x,y; X)h(y)dS g (y), D x f(x)= / d v{y) G(x,y;X)f(y)dS g (y), 

where G(x,y; A) is Green's function, 

(-A g -\)G(x,y;\) = 5 y (x). 
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Note that the single-layer potential is well-defined on and continuous across E. 
Also, D\f can be continuously extended from M^_ n * to M± and, for iGE, 



Similar, the normal derivative of S\h is continuously extended, from the left and 
right, to x E E, 



where the integrals in (2.5)-(2.6) are weakly singular ones. Also integrals in 
(2.5)-(2.6) are mutually adjoint on E in the sense of duality 



(For these results see e.g. [8]. Note that, due to the local nature of constructions 
in [8], they are valid not only for the Euclidian case considered in there but also 
for manifolds.) 

2.2. Transmission problem in time domain. In this subsection we consider 
the transmission problem for the wave equation associated with the Laplace op- 
erator, — Ag, 



where h, f e C^°(E x R). This space consists of C°°— smooth functions equal to 
for sufficiently large negative t, i.e. 

h — for t < t h , f = tort<tf. 

The results obtained below will be instrumental, although for the special case 
/ = 0, in section 4. However, as considerations for / = and / ^ are parallel, 
we will consider the general case. 

In lemma below, we use the following spectral decomposition 






(2.7) 




(2.8) 




(2.9) 
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Lemma 2.1. (Blagovestchenskii identity) Given the Cauchy spectral data (1.2) 
it is possible to evaluate the Fourier coefficients u°- ,h {t), Uj'°(t) of the waves 
u°' h (x,t), u^°{x,t), namely, 

u 0A {t)= ^ Sm ^^L~ S)) h 3 (s)ds, hj(t):= [ h^t^dSg, (2.10) 
Jt h a/ Aj Jt, 

, n /"* sin(v/A7(t - s)) f 

uf°(t) = - / lV 3 ±- ^mds, fj(t) := / fi-^d^dS,, 

Jt, v A, 



and, for \\ = 0, we should substitute sm( V s ^ fry ^ _ s y Moreover, 

w°' h (x, t) G C~(R; H\M)), u f '°(x, t) G C~(R; L 2 (M)). 
Proof. We consider first We have 

= (^(z,o,<^)) = / u°' h (x,t)<j )j (x)dv g + [ u°*( Xit ) </ ,.( x ) dVg . 

V ; JM+ J M- 

Thus, (2.8) implies that, 

d 2 t uf l {t)= [ A g u°> h {x,t)<t> j {x)dV g + [ A g u°> h (x,t)(j) J (x)dV g = 

J M+ JM- 

J (d v u^ h - d v u_l h ^ \T,4>j\T,dS g - J (u+ h - u°l h ^j \T,d v <t)j\T,dSg 

+ / u°' h (x,t)A g( p j (x)dV g = hj(t) - \juf h (t), 
Jm 

where we use that (d„u+ h — d u u°l h ^j |s = h, (u+ h — u°l h ^j |s = 0. Solving this 
second order ordinary differential equation together with the initial conditions 
u°j ,h {th) = 0, d t u°j ,h (th) = 0, provides the first formula in (2.10). Similar consid- 
erations provide the second formula in (2.10). 

To prove the second part of Lemma, we rewrite problem (2.8) in the form, cf. 
(2.3), 

%u°' h - A g u°' h = hfa, u°' h \ t<th =0. 
As, for any t, h(-, t)5 s G H^ 1 (M), we have, for any k — 0, 1, . . . , that 

oo 

X> J + l)-m(*)l 2 eC'rW- (2.H) 
J'=l 

On the other hand, representation (2.10) yields that, for j > 2, 



H h it)= f~ hS ^^d? hj (t-s)ds. 
Jo v A? 
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This implies that 

with obvious modification for j — 1. Thus, 

oo oo 

+ i)i^«? fc (or < cX> + irmwi 2 

/f oo 
^(A. + i)- 1 !^,^)! 2 ^. 
... j= i 

This inequality, together with (2.11), implies the desired result for u°' h (x,t). 
Similarly, we prove that u f >°(x, t) G C~(R; L 2 (M)). □ 

Lemma 2.1 immediately implies the following result: 

Corollary 2.2. The Dirichlet spectral data (1.3) makes it possible to find the 
trace on E x R, i.e. u°' h \j; xR , for any h G C^°(E x R). 

Proof. The result follows from formula (2.10) taking into the account that, due to 
u o,h e c™(R; H l (M)), the series (2.9) converges, for any t, in H l (M). Therefore, 
the trace 

oo 

u >%Mx,t) = J2 u °i h ( t )<l ) Mx)> 

3=1 

where the right-hand side converges in C^°(R; ff 1 / 2 (S)) □ 

Let us now warm up by considering a simple inverse problem when the part M_ 
of the manifold M is known and show how Lemma 2.1 and Corollary 2.2 can be 
used to recover M+. 

Example 2.1. Assume that, in addition to (1.3), we know the manifold (M_, g\jj_) 
Then the manifold (M + , g\M + ) can be recovered up to an isometry. 

Indeed, let h G C^°(E x R + ), i.e. h G C^°(E x R) and t h > 0. Using Corollary 
2.2, we can find -u ' h | SxR+ . 

Consider now the Dirichlet initial-boundary value problem in M_ x R + with 
Dirichlet data being - u°' /l |sxiR + - As (M_, <?|m_) is known, we can then find 
u°' h \m^xr + and, therefore, the normal derivative ctw^lsxM- Thus, we can find 

dO, hi o 0,h\ . i 

uU+ \T.xR = Oi>Uj. IsxR + Zi- 

It is shown in the proof of Lemma 4.4 below that, when h runs over C^°(E x R+), 
then M°' h |sxiR + also runs over the whole C^°(E x R+). Therefore, the set of pairs, 

{(u % xR ,d u u°+% xR ) : h G C~(E x R)}, 
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defines the graph of the non-stationary Dirichlet-to-Neumann map for M + . 

It is, however, shown in e.g. [17] that this map determines (M + , g\ M+ ) upto an 
isometry. 

2.3. Metric on E. We complete this section showing that the Dirichlet spectral 
data determines the metric on E. This result will be needed later, in sections 3 
and 4. 

Lemma 2.3. The Dirichlet spectral data (1.3) determine the distance function 
on E and, therefore, the inherited metric tensor, g\^, on E. 

Proof. Let H(x,y;t) be the heat kernel, 

d t H(x,y;t)-A g H(x,y;t)=5 y (x)5(t), x,y G M; H\ t=0 = 0. 
Then, uniformly in M x M, 

t\ogH(x,y;t) -> ~d M (x,y), 
as t — > 0, see e.g. [35]. On the other hand, for t > 0, 

oo 

H(x,y;t) = £V A 'Vj(aO0j(y), 
j'=i 

where the convergence takes place in C°°(M x M x (0, oo)). Thus, we can de- 
termine the distance dM{x,y) between any points x,y G E using the Dirichlet 
spectral data (1.3). Then the distance along E is given by 

N-l 

dx(x, y) = hminf ^ d M(zj, z j+1 ), 
j=o 

where the infimum is taken over all finite sequences of points z , zi, . . . , zn G E 
such that z = x, zn = y and d,M(zj, z j+i) < e - Having at hand ds(x,y), we can 
determine the inherited metric tensor, g\-£ on E. □ 

3. Inverse problem with Cauchy spectral data. Theorem 1.1. 

In this section we develop a procedure to reconstruct the Riemannian manifold 
(M,g) from the Cauchy spectral data (1.2). 

3.1. Prom Cauchy spectral data to the response operator. Let us consider 
the transmission problem (2.2) with A ^ o~(—A g ). Denoting, as in section 2, its 
solution by u^' h (x, A), we define the response operator by setting 

i2 A (/,/i)=«f(A)| E . 
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Lemma 3.1. The Cauchy spectral data (1.2) determine the response operator R\ 
for all\<£a(-A g ). 

Proof. Note that Green's function has the following spectral representation 

G {x ,y;X) = ±^M, (3.1) 

j=l 3 

where the sum converges the sense of operators in L 2 (M) and we assume, without 
loss of generality, that the eigenfunctions cf>j are real-valued. Thus, given the 
Cauchy spectral data and using (3.1), one can formally construct the single- and 
double-layer potentials on E, 

^(* fo ,= Jjp-^^Hym (3.2) 



, or = , / IE Mx ~' {y, ]m^ (3.3) 

• - j=l 3 



Since 



Rx(f,h) = S x h-(± + D$f, 



it looks that data (1.2) directly determines R\. However, we face the difficulty 
that series (3.1) does not converge pointwise. To deal with this difficulty, consider 
first the case when / = 0. To determine the coefficients in the Fourier expansion 

oo 

u°' h (x, A) = 5>°' fc (A), Mvwhix), (3.4) 

3=1 

we use Green's formula to get 

(«°' fc , ^)lHm) = - / + / )(A p u°' fc fc - u 0,h Ag4>j)dV g = (3.5) 

A ~ A j J M+ J M- 

- , 1 , / ~ - u+ h d u (fij\z + u°l h d u (j)j\T:)dSg = - 1 / hcfrjlzdSg. 

A ~ A j JY, a ~ A j JT, 

Since the series (3.4), (3.5) converges in fT^M), so that the trace is given by 



Wis 



u°' h (x,\)\z = -J2[ Mv)Hy)ds g (y)] 



A - A,- ' 



where the series converges in H^ 2 (E). Hence (3.2) is well-defined. To compute 
(3.2) we also need to know the Riemannian volume dS g (x) of E. By Lemma 2.3, 
it can be found from data (1.2). 
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Let us now show that (3.3) is well-defined. First note that u°' h \ M+ £ H 1 (M + ) n 
Hf oc (M + ). Thus, we can define, for any e > 0, the normal derivative d v u 0,h \^e, 
where E £ = {ieM + : d(x, E) = e}. Clearly, for any ^ G H 1 (M + ), 

/ (<9^ 'V)^^ = / A g u°> h VdV g + [ (V g u°' h ,V g V) g dV g 
= -x[ u°> h *dV g + f (V g u°>\V g V) g dV g (3.6) 

J Ml J M e + 

where = {x G M + : <i(:r, E) > e} and we have used equation (2.2). As the 
right-hand side of (3.6) has a limit, when e — > 0, and runs over if x / 2 (E) 
when \l/ runs over H\M+), this defines <9^| s . As we can choose, for any 
*|s G # 1/2 (S) its extension ^ so that < C||\I>| E ||#i/2, 

R« 0,ft (A) + | E || H -i /2(E) < C\\u > h (X)\ M+ \\m(M + ). (3.7) 

Define the sources-to-Dirichlet operator by setting 

J x h = d v u°' h (\)+\z, 
see (2.6). Taking A-derivative of (3.5), we get 



^-■u 0,/l (x, A) = £ [JjMHy)dS g (y)}^ 



n,U) (-A g - A)- V> fe (A) 



A) 2 

that converges in H 3 (M). So we have a well-defined object 

|-^(A)| S = f; [ / My)h(y)ds g (y)} d ^ {x) ^ 



d\ u v v " ^l^^v^-v^-^j (A ._ A)2 , 

where the convergence holds in if 3 / 2 (E). As -j^(J\h) = ^d u u°' h (X)\-£, we can 
compute, for any h G C°°(E), using the Cauchy spectral data of E. 

Let A G R, A 7^ Xj, and let 77 C C be the line segment from A to %T . As 

f x d 

J\h = / —{J T h) dr + Jirh, 



IT 



dr 



we have 



J\h= lim ( / —(J T h)dT + Ji T h). 



By Lemma 3.2 below, we get lim^oo Jirh = lim^oo d u u°' h (iT) + \^ = \h. This 
implies that 

-x 



1 [ A d 
J\h = -h+ lim / —(J T h)dT, 

2 J iT dr 



where the right-hand side can be computed using the Cauchy spectral data. 
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To complete the proof, we recall, see equation (2.7), that J\ is adjoint of D%. 
Thus we can find D° x f using the Cauchy spectral data. □ 

In the proof of the above Lemma we used the following asymptotics, with respect 
to singularity, of u°' h near E. 

Lemma 3.2. Let \ E C, ImX > and Re y/— A < 0. Then uniformly for 
|arg(A)| >5, 

u°' h (\)±\z^0, d u u°> h (\) ± \z ^ ±^h, as A^oo, (3.8) 
in i/ 3 / 2 (£), if 1/,2 (E) ; correspondingly. 

Proof. Let us first fix local coordinates on M near E, x — (x',x m ), where x' = 
(x 1 , . . . ,x m_1 ) are some local coordinates on E and x m is the signed distance to 
S, 

x m = ±dist(x, E) for x E M±. 
In these coordinates, we introduce 

[^le xm ^C(x m ), xEM + , 

where ((x m ) is a smooth cut-off function equal to 1 near x m = supported 
in (—a, a) with sufficiently small a > 0. Outside the a-neighborhood of E, the 
functions v± are defined to be zero. Writing 

Ag = Dim + p(x)O x m + Q(X, d x l,..., O x m-l) 

in the above coordinates and using the fact that 

f 1 \e xmV ^C(x m )\ 2 dx m < C(V^X)-\ 
Jo 

we see that (2.3), (3.9) yield 

u°' h (x, A) = v°' h (x, A) + w°' h (x, A), 
where w°' h E H 2 (M) satisfies 

(-A, - A)^ = H h (\), \\H h (\)\\ LHM) < C h {\ + IAD" 1 / 4 . 

As ||(— A g — A)" 1 )! < dist(A, a(— A g )), where the norm is the operator norm in 
L 2 (M), this implies 

\\w°' h (X)\\ H 2 {M) <C h , s \X\- 1/4 , 

when |arg(A)| > 5 > and |A| > 1. Combining this estimate with (3.9), we see 
(3.8). 

□ 
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Remark 3.1. Analyzing the behaviour of u^ h near S x R, we can show, cf. con- 
siderations leading to (3.8) and (3.9), that 

u f ' h (x,t)\ M±xR G C~(M ± x R), (3.10) 

meaning that u^' h (x,t)\ M int xR may be continued to M± x R to satisfy (3.10). 

3.2. Reconstruction of the manifold using the response operator. Recall, 
see e.g. [17], that, if (N,g), dN ^ 0, then its Neumann-to-Dirichlet operators, 
A\(N) are defined as 

A A (A0V = ^(A)U, 
where w^(x, A) is the solution to the Neumann problem 

-A g w^(x, A) = \w^(x, A), x G N mt , d u w^{x, X)\ dN = ^, 
for A ^ o"(— A^), a(-A^) being the spectrum of the Neumann Laplacian on N. 

Lemma 3.3. Given the Cauchy spectral data (1.2) it is possible to find the 
Neumann-to-Dirichlet operators A\(M±) for A ^ er(— A^) 7 where — A± stands 
for the Neumann Laplacian on M± . 

Proof. We start with A A (M_), assuming A £ (cr(-A g ) U cr(-A N ) U tr(-A^)) , 
where — A^ is the Dirichlet Laplacian in M±, correspondingly. Then, for any 
h G C°°(E), there is a unique solution, w'l(x, A) G C°°(M_), satisfying 

-A g w h _{x, A) = \w h _(x, A) in M_, d v w h _{; A)| s = -/i, 

where, as in equation (2.1), v is the unit normal pointing towards M_. 

Consider 

in M+, 
wt(x, A) in M_. 

Clearly, w h (x,\) solves (2.2) with [u^] := / = A)| E , [d^] = /i. More- 

over, with this / and /i, 

Rx(f,h) = 0. (3.11) 

These considerations show that, for any h, there is / such that (3.11) is satisfied 
and we can consider (3.11) as an equation for / when h is given. Let us show 
that the solution to (3.11) is unique if A £ (a(-A N ) U tr(-A^)). This will allow 
us to uniquely define / = f h (X) = —wt(-, A)| s as the solution to (3.11). Then, 

A x (M_)h = -f h (\). 

To prove uniqueness, assume that there is / such that 

Rx(f,0)=0. 

As A ^ cr(-A^), this implies that u^°(x, A) = 0. As [d u ii f > ] = 0, we see that 

<9^'°(-,A)| s = 0. 



w h (x,X) 
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However, A ^ a(-A N ), so that u L °(x, A) = 0, i.e. / = 0. 

Combining with Lemma 3.1, we see that the Cauchy spectral data (1.2) deter- 
mine A A (M_) for A (<t(-A s ) U v(-A n ) U tr(-A^)). Since A A (M_) is a mero- 
morphic operator- valued function with simple poles at a(-A^), this determines 
A A (M_) uniquely. 

As u_ h \s — R\(f, h) — f, we can repeat the previous arguments for A A (M + ). □ 

Theorem 1.1 follows from Lemma 3.3 taking into account that A A (M-|-) determine 
(M±, g±) up to an isometry, see [17], section 4.1. Thus to recover (M, g) we should 
just glue (M_,g_) and (M + ,g + ) along given S. 

4. Inverse problem with Dirichlet spectral data. Theorem 1.2. 

In this section, we will develop a procedure to reconstruct the Riemannian man- 
ifold (M,g) from the Dirichlet spectral data (1.3). We will assume that S C M 
consists of two open subsets Si,S 2 , Si fl Si = 0, S = S x U S 2 . As in section 
3, we assume that £ := <9S = Si U S 2 , X« = <9Sj,i = 1,2, are smooth. More- 
over, we assume that the spectra cr(— A D (Sj)), <r(— A D (M \ Si)),i = 1,2, and 
(j(-A D (M \ S)) are all disjoint. 

4.1. An approximate controllability result. Consider the following trans- 
mission problem 

{d 2 t - A g )u = 0, in (M \ S) x R, 

[u] E = 0, [d„u] E = /i G C~(E x K), (4.1) 

and denote by u(x, t) = u°' h (x,t) its solution. Note that problem (4.1) coincides 
with problem (2.8) with / = 0. 

By Lemma 2.1, u°' h G C^°(M; H l (M)) and we can define an operator 

W : C^°(£ x R) — > if 1 (M), W7i := u°> h (0), 
which is called the wave operator associated with problem (4.1). 
Theorem 4.1. Let a(-A D (S)) n a(-A D (M \ S)) = 0. T/jen toe set 

Y = {Wh : /i G C~(£ x R)} (4.2) 

is dense in H l (M). 

Proof. Assume that ip G {^(M))' = H' 1 (M) is orthogonal to Y, 

(u^ h (0)^) H i (M)xH - HM) = (4.3) 
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for all h G C+(E x R). Let e be the solution to the problem, 

e tt - A g e = 0, in M x R, (4.4) 
e|t=o = 0, e t \ t=0 = ip. 

Then similar considerations to those at the end of proof of Lemma 2.1 show that 
the weak solution, e(x, t) of (4.4) satisfies, 



e(x,t) =2_^ 7= (4>,(j>j)H-i(M)xHi(M)(i>j +t{^Al)H-^M)y.m(M)<Pl (4.5) 



3=2 



and e G C(R;L 2 (M)) n C\R; H'^M)). Observe that, as < C(\j + 

l) 1 / 2 and Aj > Cj 2 l m , e has a well-defined trace, in x R), on E x R with 



J (^,0j)0j|s + -» e(g,t)| SxR , as J -> oo, (4.6) 



j=2 V 

in x R). 

Let us show that e(x, t)|sxR = 0. Choosing /i G Co°(E x R_) and using Green's 
formula, we obtain from (4.1)-(4.4) that 

</MxIR_ 

= {u ' h (T),ilj) H i iM)xH -i (M) + hedS g dt = hedS g dt. 

This yields that e| Sx (_ OOi0 ) = 0. As by (4.4) e(x, s) = —e(x, — s), we see that 

supp(e| SxM ) = E x {0}. (4.7) 

Next we show that 

e|Ex(-i,i) G ^ 1/2 (S x (-1, 1)) := (H 1 ' 2 ^ x (-1, 1)))' . (4.8) 

Let X be a local, near S, vector field on S 1 such that X\% = d v \z. Let h G 
# 1/2 (£ x (-1, 1)) and H G H l (S x (-1, 1)) be its continuation into S x (-1, 1), 
such that H = outside the domain of definition of X and 

||#||ffl < C||/l|| H l/2. 

Denote by E the primitive, with respect to t, of e in 5 x (—1, 1), 



( /X-t') t 2 
= ^. * ' (MM*) + o WMMx), (4-9) 

j=2 J 



INVERSE SPECTRAL PROBLEMS ON A CLOSED MANIFOLD 15 
E E C(R; H^S)) n C\R; L 2 {S)). Integrating by parts, we get 



/ hedS g dt= I (H ■ Xe — X C H ■ e) dV g dt (4.10) 

isx(-l,l) JSx(-l,l) 

= [ ([H-XE}\ t=l -[H-XE}\ t= ^)dV g - [ (d t H-XE + X c H-e)dV g dt, 

JS JSx(-lA) 



where X c is the first-order operator adjoint to X. By (4.5), (4.9), the right-hand 
side of (4.10) can be estimated by 

C\\H\\ H i( Sx (-i,i)) < C / ||^||i/i/2(sx(-i,i))- 

Thus the left-hand side of (4.10) is bounded for any h G H 1 ' 2 ^ x (-1,1)), 
proving (4.8). 

Now (4.7) implies that e(x, £)|exm = J2i=o ^i{x)d l 8{t), with some finite /, see e.g. 
[12, ex. 5.1.2]. Thus, (4.8) yields 

eOM)|s X R = 0. (4.11) 

The last step of the proof is to show that this equation yields that e = in M x R. 
Using relation (4.6), equation (4.11) and making the partial Fourier transform, 
t — > k, we see that the distribution, e(x, k) e «S'(E x R), satisfies 

oo 



,l^(S(l;-i/Ai) 
e(z,*0 = HoL 7?= WWjls 

2 i=2 V A i 

1 (^f /c I /X - ) 

+^(A:)(^,0 1 )0 1 |e--^ 7f - (^0 J -)0 J -|s)=O- 

This implies that (■?/>, <Pi)h- 1 (m)xh 1 (m) = and, for any j = 2, . . . , 

X] (^^i)H-i(M)x^i(M)0ils = 0, (4.12) 

3 'J 

where the last sum takes into account eigenspaces of an arbitrary multiplicity. 
Consider the function 

$(z) := ^ (^,0j)H-i(M)xHi(M)0i(^), x E M. 

3 3 

It satisfies the Dirichlet boundary condition, $|s = 0, see (4.12), and, as 4>j are 
eigenfunctions of — A g with Xj = Aj, the equation 

-A g $(z) = Aj$(x), x G M. 

Thus $|s is an eigenfunction of — A D (S), while $|m\s is an eigenfunction of 
-A D (M\S). However, as a(-A D (S)) Da(-A D (M\S)) = 0, we have that % = 
or $|m\s = 0. In any case, by the uniqueness of zero-continuation for elliptic 
equations, this yields that $ = everywhere in M. As different (f>j, corresponding 
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to Xj = Ay, are linearly independent, this implies that (ijj, (f)j)H- 1 (M)xH 1 (M) = ^ 
for all j = 1, 2, . . . . Thus, e = in M x R and, therefore, ^ = 0. □ 

4.2. Approximate controllability with given trace at final time. In this 
section we denote E to be either Ej, % = 1,2, orE. Lemma 2.1 makes it possible 
to introduce a quasinorm 

oo 

\h\ 2 := \\Wh\\l 1{M) = £(A, + 1)|^(0)| 2 . (4.13) 

j=i 

It is classical for the control theory, see e.g. [20] or [22] in the context of inverse 
problems, to define the space D 1 of the generalized sources by introducing the 
equivalence relation, 

h= E h if u°' h (0) = u°' % (0), 

and completing C^°(E xR)/E with respect to (4.13), 

D 1 : = cl (C7(E x R)/£) . 

Then, by Theorem 4.1, we can extend the wave operator W, see (4.2), from 
C~(E x R) onto D\ 

Wh:=u°' h (0), W:D 1 ^H\M), 
unitary operator. 

Moreover, as Wh = X^=i rf^j -f^(M) and the Fourier coefficients k^, for any 
h G D 1 , can be explicitly evaluated using the Dirichlet spectral data (1.3), see 
the first formula in (2.10), it is possible to find, for such ft,, the trace, 

oo 
3=1 

The above considerations give rise to the following lemma. 

Lemma 4.2. Assume a(-A D (S)) D a(-A D (M \ S)) = 0. Then the subspaces, 
Dl : = {feGi) 1 : Wh\^ = 0} C D 1 , 

are uniquely determined in terms of the Dirichlet spectral data (1.3). 
Moreover, the wave operator W , restricted to D~, 

W-.D^^H 1 -, H±:={aEH l (M):a\z = 0}, 

is unitary. 
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4.3. Finding eigenvalues and eigenfunctions in subdomains. In this sub- 
section we denote by S one of the manifolds Si, % — 1,2, M \ S i: % — 1,2, M \ S 
and by X n {S), <f) n (-;S) we denote the eigenvalues and ortho normal eigenfunctions 
of — A D (S). By the max-min principle, 

X n (S) = max min(V 9 M n , V 9 M n ) L2( ^, 

Ul,...,U n _l u n v I 

where the maximum is taken over u\, U2, ■ ■ ■ , u n -\ e H\{$) and the minimum is 
taken over w n G Hq(S) that satisfies 

(Mn,M P ) L 2 ( s) = 0, p= l,2,...,n- 1, (tip,Up) £ a ( g) = 1, p=l,2,...,ra. 

The minimizer tt n (:r) is then an normalized eigenfunction corresponding to the 
eigenvalue X n {S). Now consider the following max-min problem 

oo 

max min(V 9 M„, V S M„) L 2 (M) = max min V" Xj \u njj \ 2 , (4.14) 

U1,...,U„_1 U n Ul,...,U n -l Un * ' 

where u n j are the Fourier coefficients of u n , i.e. u n (x) = YLjLi u n,j<f>j( x )- Here 
the maximum is taken over u±, u-i, ■ ■ ■ , w n -i € H X [M) satisfying u p \^ = and the 
minimum is taken over u n G if 1 (M) satisfuing u n |g = and 

oo 

(u n ,u p ) L 2 {M) = ^UnjUpj = 0, p = 1,2, . . . ,n - 1, (4.15) 
i=i 

oo 

(up, u p ) L 2 {M) = ^2 \ u P,j\ 2 = !> p = 1, 2, . . . , n. 
j'=i 

Then t„(E) are the eigenvalues of the Dirichlet Laplacian on the direct sum of 
L 2 (S) and L 2 (M\S), so that 

WE)}^ = a(-A D (S)) U ( r(-A D (M \ 5)). 

The sequence of the corresponding minimizers, u n (x; S) consists of orthonormal 
eigenfunctions of this operator. However, due to the assumption er(— A D (S)) fl 
cr(— A D (M \ 5)) = 0, any such eigenfunction is equal to on M \ S or S. Thus, 
any w n (x; E) is either an eigenfunction of — A D (S) extended by to M\ S, or an 
eigenfunction of —A D (M \ S) extended by to S. 

On the other hand, Lemma 4.2 together with equation (2.10) make it possible to 
evaluate the right-hand sides in (4.14), (4.15) using the Dirichlet spectral data. 
This leads to the following result 

Lemma 4.3. Let a(-A D (S~)), a(—A D (M \ $)), i = 1, 2 and a(-A D (M \ S)) 
be all disjoint. Then the Dirichlet spectral data (1.3) determine uniquely the 
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eigenvalues X n (S), n = 1,2, ... , S = S, Si, M \ Si, % = 1,2, and M \ S. They 
determine also the the generalized sources h n (S) such that 

Wh n (S) = u°> h ^\x,0) = \ MX; ^' f ~ S> (4.16) 
y ! V ' ; (0, x G M\S. V ; 

In addition, the Dirichlet spectral data determine the Fourier coefficients of the 
extended eigenfunctions Wh n (S)(x), 

oo 

Wh n (S) = "nAS)M x l x E M. (4.17) 
3=1 

We note that these sources h n (S) are determined up to a unitary transformation 
in the eigenspace corresponding to the eigenvalue X n (S). 

Proof. Recall that, by formula (2.10) we can evaluate the Fourier coefficients 
u®' (0) for any h G D 1 . Thus, we can evaluate 



H n {T) = max min V A>J' l "(0)| 2 , 

h\,...,h„-\ h n — 



3=1 

where the maximum is taken over hi, h 2 , . . . , h n -\ G D~ and the minimum is 
taken over h n G D~ with 

oo oo 

j>j*»(0)tij**(0) = 0, p = 1, . . . , n - 1, \uT P ( )\ 2 = 1, P = 1, • • • , n. 

j=l j=l 

It follows from Lemma 4.2 that, for any n — 1,2, ... , /i n (S) = t n (E), providing 
a(-A D (S)) n a(-A D (M \ S)) = 0. 

Repeating this construction with S equal to S, Si, and S 2 , we obtain the sets 
a(-A D (S)) U a(-A D (M \ S)) = a(-A D (S 1 )) U a(-A D (S 2 )) U a(—A D (M \ S)), 
and 

a(-A D (S 1 )) U a(-A D (M \ S 1 )), a(-A D (S 2 )) U a(-A D (M \ S 2 )). 

Asa(-A D (S~)), a(-A D (M\Si)),i = 1, 2 and a(-A D (M\S)) are all disjoint, by 
intersecting the above sets we find the desired eigenvalues X n (Si), n — 1, 2, . . . , % — 
1, 2 and X n (M \ Si), n = 1,2, ... ,i = 1,2, as well as X n (M \ S), n = 1, 2, . . . . 

Then, identifying the corresponding subsequence of t n (S) and related generalized 
sources /i n (S), we determine, for each S, the generalized sources h n (S) such that 
Wh n (S) are equal to the extended eigenfunctions (4.16). 

Recalling formula (2.10), we prove the last part of the Lemma. □ 
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4.4. Inverse problems in subdomains. Proof of Theorem 1.2. Our proof 
of Theorem 1.2 is based on Lemma 4.3. Namely, we will show that, having 
at hand the eigenvalues X n (Si), i — 1,2, and A n (M \ S) and also the Fourier 
coefficients, K n j(Si), i — 1,2, K n j(M \ S), it is possible to determine, up to an 
isometry, the Riemannian manifolds (Si,g), i — 1,2, (M \ S,g). Gluing them 
along Ej we recover (M,g). 

Recall that if, for (N,g), dN ^ 0, we do know its Dirichlet eigenvalues X n (N) 
and traces on ON of the normal derivatives of the eigenf unctions, d u ip n \sN, then 
these data determine (N,g) up to an isometry, see e.g. [17]. However, in the 
case of the Dirichlet spectral data, we have only the Dirichlet values 4>j\n an d, 
moreover, the convergence of the Fourier series (4.17) is only in H 1 (M) preventing 
us from identifying d u <f> n (S)\^. Therefore, we will use another approach within 
the BC-method, described in section 4.1 of [17]. To explain it, consider the 
intial-boundary value problem in N x R, 

(d?-A g )w F = 0, inivxR + , 

W F \ 9 NxR + = F, w F \ t=0 = 0, d t w F \ t=0 = 0, 

with F G C+(dN x R + ). The energy, at time t, of the wave w F is then defined 

as 

E(w F ,t) =\ [ (\d t w F (x,t)\ 2 + \V g w F (x,t)\ 2 g )dV g (x). 
1 Jn 

It is shown in [17] that, given the energy flux 

11(F) := lim E(w F ,t), 

for any F G C^°(dN x R+), it is possible to determine {N,g) up to an isometry. 

Therefore, Theorem 1.2 is an immediate corollary of the following Lemma, 

Lemma 4.4. Let E = E 1 UE 2 divide M into regions Si, i = 1,2, andM\S. Con- 
sider the initial-boundary value problems (4-18) with N equal to Si, S% and M\S. 
Then, assuming that a(-A D (S~)), a(-A D (M\Si)),i= 1,2, and a(-A D (M\S)) 
are all disjoint, the Dirichlet spectral data (1.3) uniquely determine the energy 
flux 11(F) in each of these subdomains. 

Proof. We start with the case iV = S\. Then (4.18) takes the form 

(d*-A g )w F = 0, inSixR + , 

w f \j: iX r + = F, w F \ t=0 = 0, d t w F \ t=0 = 0, 

with F G C^°(E 1 x R+). Let us first show that, for any such F, there exists a 
unique h F G C^°(E 1 x R + ) such that w F = u°' h \ Sl , where 

h=[ hF > (4.20) 
0, x G E 2 . 
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To this end we consider, in addition to (4.19), the problem 

{d 2 t - A g )w c F = 0, in (M \ S ± ) x R+, 
w c f \ SiX r + = F, w c F \ t=0 = 0, d t w c F \ t=0 = 0, 
and introduce the function u, 



u= ,w F , xeS,xR + , 

^w c F , xe(M\ S 1 ) xR + . K J 



Then u solves the transmission problem (2.8) with / = and 
h = 



h F = (d u w c F - d v w F ) |eixr+, OM) G Si x R + 
0, (x,t)GS 2 xM+, 

i.e. u = u 0,h . Using considerations similar to those in the proof of Lemma 3.3, 
we show the uniqueness of such h. 

By Corollary 2.2, we can then find, for h G C^(Iq x R+), 

Ai/i = u°> h | ElxR+ . (4.22) 

As shown earlier, for h = on E 2 , the operator Ai is an invertible operator 
in C^°(Si x R+). Thus, we can use equation (4.22) with the right-hand side 
being F, to uniquely determine h F . Observe that the extended eigenfunctions 
{W (/i n (5i))}~ =1 together with {W {h k {M \ Si))}^ form an orthonormal basis 
in L 2 (M). Thus, using (4.17), we can evaluate the Fourier coefficients of u°' h (-, t) 
with respect to this basis, 



u° 



> h (x, t)=J2 W (M&)) ^ + E W (h k (M \S 1 ))(x), 



n=l k=l 



where the index F indicates that h is of form (4.20). This expansion, together 
with the definition of u, see (4.21), shows that 

oo 

w F (x, t) = e&Wix, = E w nAt) w (h n (Si)) (X), 

n=l 

oo 

\Mt)\\h iSl) = Y,\ w nAt)\ 2 , (4-23) 

n=l 

where 6i is the characteristic function of Si and we identify functions in Si with 
their extensions by to M \ Si. 

Similarly, 

oo oo 

d t W F (x,t) = J2 d t W n, F (t) W (h n (Si)) (X), 11^(011^(50 = E \ d t W n,F(t)\ 2 
n=l n=l 
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The last step of the construction is based on the observation that, when F G 
Cq°(£i x R + ), then w F G Hq(Si) for large t and, therefore, the eigenfunction 
expansions (4.23) converge in H l (S\). Thus, we can find 

oo 

l|v s M-,*)||i 2(5l) = USi)\wnAt)\ 2 - 

n=l 

Combining the above two equations, we find the energy flux 11(F). 

Similar considerations, with S2 and M \ S2 and M \ S and S instead of S± 
and M \Si, show the possibility to evaluate the energy flux when N = S2 and 
M\S. □ 



5. Some remarks on further generalizations and open problems 

Remark 5.1. The inverse problem with the Cauchy or Dirichlet spectral data on a 
compact Riemannian manifold without boundary studied in Sections 3, 4 can be 
generalized to the problem when M has boundary dM 7^ with e.g. Neumann 
(or Dirichlet) boundary condition. Then (1.2) and (1.3) consists of the Neumann 
eigenvalues of M and the Cauchy or Dirichlet data of the eigenfuctions on a 
closed hypersurface £ C M, dM n £ = 0. The methods to solve these problems 
are very similar to the described and are based on the version of the BC-method 
with data on a part of boundary, see section 4.4 [17]. 

Remark 5.2. Consider the non- stationary Green function, G(x, y; t), for the acous- 
tic wave equation in (M,g), where M may have a non-trivial boundary, 

(d 2 t - A g )G(x, y; t) = S y (x)S{t), in M x R, 
G(x, y; t)\ dM xR+ = 0, G(x, y; t) | t<0 = 0, 

where y E £, £ C M such that £ fl dM = and the boundary condition is 
void when dM = 0. It often happens in practice, for example in geophysics or 
ultrasound imaging, that one can measure G(x,y;t) for x again running only 
over £. Thus, we come to the inverse problem with non-stationary data being 
G(x,y;t), x, y G £, t > 0. 

Taking the Fourier transform of the given G(x, y; t) in t, we find the Green func- 
tion G(x, y; k), cf. [15], 

(-k 2 - A g )G(x,y; k) = 5 y (x), G(x,y;t)\ 9M = 0. 

Thus, from practical measurements we can find G(x,y; k) for x, y G £, k G R. 

Note that G(x,y;k), x,y G £, is the integral kernel of a meromorphic, with 
respect to k G C, operator- valued function in £ 2 (£). In terms of the eigenvalues 
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and eigenfunctions of — A g , 

-\ ^ 

C(x,y;A;)| SxS = ^ — _ - ^-(a;)| s 3 -(y)| s , 

where the convergence in the right-hand side is understood in the sense of oper- 
ators in £ 2 (£). So finding the poles of G(x,y; k) is equivalent to determination 
of y/Xj. At the pole \fXj, the residue is given by the integral operator with the 
kernel, 

Ii£s(G(x,y;-),y/Xj) = -^= ^ M x )\vMy)\v- 

The knowledge of this kernel allows us to find the functions (f>j\z up to an orthog- 
onal transformation in the eigenspace corresponding to Xj, see [15], [31]. 

Therefore, the dynamic inverse data G(x,y;t), x, y G X, t > makes it possible 
to find the Dirichlet spectral data on E. 

Remark 5.3. As shown in section 4, the Dirichlet data (1.3) determine (M,g) 
when the spectra a(-A D (Sl)), a(-A D (M \ S l )),i = 1,2, and a(-A D (M \ S)) 
are all disjoint. It is interesting to understand whether this condition can be 
removed. 

In general, it is important to find if the Dirichlet spectral data (1.3) determine 
(M, g) even in the case when E is connected. 
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